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Total and parity-projected level densities of iron-region nuclei are calculated mi- 
croscopically by using Monte Carlo methods for the nuclear shell model in the 
complete (pf + 039/2 )-shell. The calculated total level density is found to be in 
good agreement with the experimental level density. The Monte Carlo calculations 
offer a significant improvement over the thermal Hartree-Fock approximation. Con- 
trary to the Fermi gas model, it is found that the level density has a significant 
parity-dependence in the neutron resonance region. The systematics of the level 
density parameters (including shell effects) in the iron region is presented. 



1 Introduction 

Neutron-capture reactions play an important role in nucleosynthesis, e.g. in 
the s and r processes. Their rates are strongly affected by the corresponding 
nuclear level densities around the neutron resonance region. Most conventional 
calculations of the jnuclear level density are based on the Fermi gas model (e.g. 
the Bethe formula)El. A phenomenological modification is often adopteda where 
the excitation energy E x in the Bethe formula is backshifted by A, giving a 
total level density of 

p{E x ) « 9 ^§a-i(E x - ArieV^M (1) 

with g = 2. The backshift A originates in pairing correlations and shell effects, 
while the parameter a is determined by the single-particle level-density at the 
Fermi energy. By adjusting the value of a for each nucleus, Eq. (0) describes 
well a large volume of experimental data. However, the Fermi gas model grossly 
underestimates the value of a, and consequently it is difficult to predict the 
level density to an accuracy better than an order of magnitude. Much less is 
known about the parity-dependence of the level density. 

We have been studying the level densities of iron-region nuclei! -by applying 
the recently proposed shell model Monte Carlo (SMMC) methoda within the 
the full pf- and Ogg M-shell. As a finite temperature method, it is particularly 



suitable for calculations of level densities (see below) . By introducing parity- 
projection methods for the auxiliary fields we are also able to calculate parity- 
projected level densities. We remark that the (pf + Ogg/2) model space is 
sufficiently large to describe the important excitations around the neutron- 
resonance energies in this mass region [E x ~ 5 — 15 MeV). 

2 Calculations of level densities 

We adopt an isoscalar Hamiltonian of the formi 

H = J2 tana + g P {0 ^ • - x £ k x OW ■ . (2) 

a X 

The single-particle energies e a are determined from a Woods-SaxorLpotential 
V plus spin-orbit interaction with the parameters quoted in Ref. El. p(°' 1 )t 
denotes the T = 1 pair-creation operator and O^-^ are multipole operators 
with radial part of dV/ dr EL In the present calculations we include multipole 
terms with A — 2,3 and 4. Core polarization effects are taken into account 
through the use of renormalization factors k\. We use &2 = 2, = 1.5 and 
£4 = 1. The pairing strength go is determined from the experimental odd-even 
mass differences for spherical nuclei in tha-mass region A = 40 — 80. This 
Hamiltonian satisfies the modified sign rule til, and therefore has a good Monte 
Carlo sign for even-even nuclei. This enables us to perform accurate Monte 
Carlo calculations. 

In SMMC the canonical thermal energy E((3) = (H)p is calculated as 
a function of inverse temperature (3. Particle-number projection (for both 
protons and neutrons) is implemented exactly!! The canonical partition func- 
tion Z(j3) is determined by a numerical integration of E(/3), In [Z {(3) / Z ({))] = 
— Jq d/3'E((3'). The level density p(E) is then calculated in the saddle-point 
approximation from 

p{E) = (2ir[3- 2 C)- 1/2 e s ; 

S(E) =(3E + \n Z(f3) , (3" 2 C{(3) = -dE/d/3 . (3) 

C is the heat capacity calculated by numerical differentiation of E({3). To 
compare with experimental data, it is necessary to find p as a function of the 
excitation energy. For that purpose we determine the ground-state energy by 
extrapolating E((3) to (3 — * 00. 

To calculate the parity-dependence o£the level density, we have introduced 
parity-projection techniques in SMMC 13. Using the Hubbard-Stratonovich 
representation for e~" H and the projection operators P± — (1 ± P)/2 (P is 



the space reflection operator), we can write the projected energies E±{0) = 
Tv(HP ± e- fJH )/Tr(P ± e- (3H ) in the form 



The integration over the auxiliary fields a is performed with the usual Monte 
Carlo weight function W(a) — G(a)((a), where G is a gaussian factor a: 
£(<r) = Tr U a is the partition function of the non-interacting propagator U ( 
Both U a and P can be represented by matrices in the single-particle space. 
E±{(3) are then calculated through matrix algebra in this single-particle space, 
similarly to the calculation of E((3). Once we find E±((3), the calculation of 
the projected densities p±(E) is analogous to that of the total level density. 



3 Results 

3.1 56 Fe 

Results for total level density in 56 Fe are presented in the left panel of Fig. |l| as a 
function of E x . The SMMC results are the solid squares and include statistical 
errors (often the errors are too small to be visible). Although it is difficult 
to measure the total level density directly, it can be reconstructed from a few 
parameters that are determined experimentally. The solid line in Fig. [I] shows 
this experimental level density as determined from charged particle spectralij. 
Our SMMC result is in good agreement with the experimental results. We 
have observed similar agreement (to within a factor of two) for other nuclei in 
the iron region. Consequently, accurate level densities can be calculated in the 
present approach. In particular, we extract the level density parameters a and 
A via a fit of our microscopically calculated level densities to Eq. ([[]). Using the 
energy range 4.5 MeV < E x < 20 MeV, we obtain a = 5.780 ± 0.055 MeV -1 
and A = 1.560 ± 0.161 MeV for 56 Fe. In Fig. | we also compare our SMMC 
results with those of the thermal Hartree-Fock approximation (HFA)El (dashed 
line). The kink observed in the HFA level density around E x ~ 9 MeV is a 
signature of a shape phase transition, but it is washed out in the SMMC due 
to strong two-body correlations. 

The Fermi gas model predicts equal positive- and negative-parity level 
densities at all energies. However, this is unrealistic in the neutron resonance 
regime of iron-region nuclei where the excitation energy is comparable to the 
energy gap among major shells (i.e. the gap between the pf and g 9 / 2 shells). 
The SMMC parity-projected level densities of 56 Fe are shown in the right 
panel of Fig. [p. They are well fitted to Eq. ([!]) with g = 1, but with parity- 
specific parameters a± and A±. We find a + = 5.611 ± 0.073 MeV -1 , A + = 
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Figure 1: Level densities of 56 Fe. Left: Total level densities in SMMC (solid 
squares) and in the HFA (dashed line). The solid line is the experimental level den- 
sityllj. Right: positive-parity (circles) and negative-parity (triangles) level densities 
in SMMC. 

0.550 ± 0.196 MeV and a_ = 6.209 ± 0.625 MeV -1 , A_ = 3.172 ± 1.637 MeV. 
Since negative-parity states in 56 Fe are possible only when g 9 / 2 is populated, 
p_ is lower than p + at low energies. Thus the backshift A_ is larger than A+. 
On the other hand, at high excitation energies positive- and negative-parity 
level densities are approximately equal, resulting in a_ > a+. 

3.2 Systematics and shell effects 

We discuss now the nucleus-dependence of the level density parameters a and A 
for even-even nuclei in the 50 < A < 70 region: 54 " 58 Fe, 58 ~ 64 Ni and 64 ~ 68 Zn. 
Among them 54 Fe and the Ni isotopes have /7/2-subshell closure for protons 
and/or neutrons (Z or N = 28). Fig. || shows the calculated values of a. It 
is interesting to see whether shell effects (at Z or N — 28) can be observed 
in the level density parameters. We have found enhancement of the backshift 
parameter A at Z or N = 28 El. However, no strong shell effects are observed 
in the single-particle level density parameter a and it increases smoothly as a 
function of A. For the parity-projected level densities we find that a_ increases 
as a function of A more moderately than a + . This is because more excitations 
to the 59/2 orbits become possible (for fixed E x ) as A increases. It would 
be interesting to investigate how the parity-dependence affects the neutron- 
capture reaction rates. 
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Figure 2: The single-particle level density parameter a as a function of A, for the total 
(left) and parity-projected (right) level densities. Points corresponding to isotopes 
with a given Z are connected by dashed lines. In the right panel, values of a+ arc 
shown by circles and values of a_ by triangles. 



4 Conclusion 

We have used the auxiliary fields Monte Carlo method to calculate the level 
density of iron-region nuclei in the complete (p/ + 0g9/ 2 )-shell, and found good 
agreement with experimental results. The introduction of a parity-projection 
technique in the SMMC allows us to study the parity-dependence of the level 
density parameters. The systematics of these parameters have been presented. 
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